- 1 - 



Os 
OS 



O 



> 
o 

oo 
O 
O 

so 
OS 

a: 
6 

CO 

c3 



13 




THE UNIVERSITY OF TOKYO 



4ggi Si? ^search Center for the Early Universe RES CEU-25/96 



UTAP-242/96 



Re-examination of the Hubble Constant 

from the Sunyaev-Zel'dovich Effect: 
Implication for Cosmological Parameters 



1,4 



Shiho Kobayashi ' , Shin Sasaki and Yasushi Suto 

1 Department of Physics, The University of Tokyo, Tokyo 113, Japan 
2 present address: The Racah Institute for Physics, 
The Hebrew University, Jerusalem 91904, Israel 
3 Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-03, Japan 
4 RESCEU (Research Center for the Early Universe), School of Science, 
The University of Tokyo, Tokyo 113, Japan 

e-mail: shiho@alf.fiz.huji.ac.il, sasaki@phys.metro-u.ac.jp, suto@phys.s.u-tokyo. ac.jp 

Abstract 

We performed a systematic and comprehensive estimate of the Hubble constant H via 
the Sunyaev-Zel'dovich temperature decrement and the X-ray measurements for clusters 
up to the redshift z = 0.541, with particular attention to its dependence on the density 
parameter Q Q , and the cosmological constant Xq. The resulting values of H are largely 
consistent with that derived from the Cepheid distance for nearby galaxies. We discuss 
how the angular diameter distance d A vs. z diagram from the accurate SZ measurements of 
high- z clusters can distinguish the values of fl and A . 
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1. Introduction 

There have been a lot of attempts to determine the Hubble constant Hq 
(= 100/i kms -1 Mpc -1 ) from the Sunyaev-Zel'dovich (SZ) temperature decrement 
and the X-ray measurements (Sunyaev & Zel'dovich 1972; see Rephaeli 1995 for a review). 
The current results are subject to several uncertainties due to simplified model assumptions 
in describing the realistic cluster gas distribution as well as the limited spatial resolution of 
the X-ray and radio observations (Birkinshaw, Hughes & Arnaud 1991; Inagaki, Suginohara 
& Suto 1995). Nevertheless this method is very important for two reasons; it does not rely 
on any local (or empirical) calibrator, and it can be used even at high redshift z ^ 1 in 
principle where no other reliable distance indicator is known, perhaps except for the SNe la 
(Perlmutter et al. 1996). 

As discussed by Silk & White (1978) and Birkinshaw et al. (1991), the observation of 
the SZ effect determines the angular diameter distance d A (z) to the redshift z of the cluster. 
For z <C 1, d A (z) is basically given only in terms of H (~ cz/H , with c being the light 
velocity). If z ^ 0.1, however, the density parameter fl , the dimensionless cosmological 
constant A , and possibly the degree of the inhomogeneities in the light path (see §3 below; 
Dyer & Roeder 1972) make significant contribution to the value of g?a(<2) as well. 

In addition, the current estimate of Hq from optical observations is converging with 
fractional uncertainty of about 10 percent level, h = 0.7 ±0.1 (Freedman et al. 1994; Tanvir 
et al. 1995). This implies that one may place useful limits on Qq and Ao if accurate and 
reliable SZ observations are carried out even for a single high-z cluster. Unfortunately the 
non-spherical distribution and dumpiness of the cluster gas among others would inevitably 
contaminate any SZ observation (Birkinshaw et al. 1991; Inagaki et al. 1995), and this is 
why we have to proceed statistically. In this Letter, we present estimates of H for seven SZ 
clusters (Table 1) by taking proper account of the effect of f2 and A for the first time. 

2. Angular Diameter Distance from the SZ Effect 

The temperature decrement in the cosmic microwave background (CMB) by the SZ 
effect observed at an angle 9 from the center of a cluster is given by the integral over the 
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line of sight: 
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where T gas and n e are the temperature and the number density of the electron gas in the 
cluster, To is the temperature of the CMB (= 2.726K) and x is the dimensionless frequency 
h p T>/kBT (y is the frequency of the radio observation, h p is the Planck constant), and the 
other symbols have their usual meanings. The relativistic correction to the above formula 
is negligible for T gas £ 15keV and x £ 1 (e.g., Rephaeli 1995). 

Similarly the X-ray surface brightness S x of the cluster observed at 9 in the X-ray 
detector bandwidth v\ ~ z/ 2 is written as 

Ssm = —^ i r dl r' 1+ "^, (3) 

4Tr(l + zp J-oo Jui{i+z) dudV 

where z is the redshift of the cluster, and d 2 L x / ' dvdV is the X-ray emissivity per unit 
frequency. If T gas is greater than several keV, the line contribution is negligible^ and S x is 
well approximated by that of the thermal bremsstrahlung: 

l^T = ^)nl Vff (T^ V )e^(--^y (4) 

(T \ ^ ( 271 2 M 

a[9as) ~ 3m e c*{3m e c*k B T gas ) l + X' ^ 

We assume the primordial abundances of hydrogen and helium with the hydrogen mass 
fraction X = 0.755. The velocity averaged Gaunt factor of the thermal bremsstrahlung, 
gff(T gas , u), is computed according to Kellogg, Baldwin & Koch (1975). 

The conventional assumption is that the cluster gas follows a spherically symmetric 
isothermal /3-model: n e (r) = n e0 /(l + r 2 j ' d\{z)9l)^' 2 , where 9 C and n e0 are the angular core 
radius (fitted from the radio or X-ray imaging) and the central density of the cluster gas, 



1 The X-ray flux due to line emission for a solar abundance plasma at T gas = 6.8keV, for 
example, is about 6 percent of that of bremsstrahlung in the ROSAT band. 
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and d A (z) denotes the angular diameter distance to the redshift z of the cluster. Then one 
can integrate equations ([!]) and (Q) over the line of sight, and their central values (9 = 0) 
reduce to 

AT(0) fc B T gas IXW2 - 1/2) 

— = -2e(^)^T mec2 vm/2) n e0 d A ( Z )9 r , (6) 
_ a(T 3a5 )fc B T gas n 2 e0 d A (z)6 x T(3/3 X - 1/2) rWH*)/^ 

xl J " 4^(1 + ^ r(3/3 x ) J hp » l{ i +z)/ k B T 3as 9ffG ax ' Uj 

where r(x) denotes the gamma function. Although the values of (3 and 6 C obtained by the 
radio (j3 r and 8 r ) and the X-ray observations (/3 X and # x ) should be identical as long as the 
spherical /3-model is exact, different values are often reported in the literature. While we 
distinguish them in the above equations following the literature, but the difference due to 
the different choices of the parameters should be regarded as an additional uncertainty in 
the present context (see §3 below). 

Finally equations (|6|) and (^) can be solved for d A : 

a(T gas )k B T gas ( m e c 2 V T(3(3 X - 1/2) / T(3(3 r /2) \ 2 
l^h p a 2 T {l + zY\k B T gas ) r(3/3 x ) \nW r /2 -1/2) ) 
(AT/T(0)) 2 6 X rV*(i+*)/*BT ff „ _ 

£ 2 (V)S X (0) 6 2 JhpUiil+zj/kBTgas gffG 

Comparison of the above estimate with the theoretical formula for d A in the Friedman- 
Robertson- Walker model, yields a relation among H Q , Q and Ao (e.g., Silk & White 1978; 
Inagaki et al. 1995) as we discuss extensively below. 



x ■ 



3. H , Q AND A from a Sample of SZ Clusters 

In order to estimate d A and see the constraints on Hq, Qq and Ao, we use seven 
clusters of which we can collect the necessary data in radio and X-rays for the analysis 
from the literature and the private communications^. In the case of CL0016+16, we adopt 
the value of 5" x (0) given in White, Silk & Henry (1981) and Birkinshaw, Gull, & Moffet 

2 We did not include A665 (Birkinshaw et al. 1991) since we could not find the proper 
X-ray data from the literature. 
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(1981). Otherwise we first compute S x (0) in equation (0) from the value of n e0 given in the 
literature of the X-ray observation, and then estimate d A (eq.[§]) so as not to include the 
line emission contribution (even though it is small). Our estimates of H Q listed in Table 1 
are for the Qq = 1 and Ao = model, just for example; they are in good agreement with 
the previous work as long as we adopt the same parameters. In what follows, we assign la 
errors simply by taking account of those of n e o, T gas and AT(0) in quadrature. The errors 
are dominated by the uncertainty in the SZ temperature decrement. In reality those values 
would be strongly dependent on the fitted parameters (3 and 6 C , and thus it is likely that 
our quoted error bars underestimate the actual ones. 

Figure [l] illustrates the dependence of the estimated H on Q for the two highest- z 
clusters in Table 1. More specifically, for CL0016+16 (z = 0.541) and A2163 (z = 0.201), 
h = 1.03±8:|| and h = 0.82±g;|| in the fi = 1 and A = model, and h = 1.43±g;|| and 
h = 0.94±g in the Qq = and Ao = 1 model, respectively. Figure || summarizes the 
estimated Hq for all clusters in models with four different sets of f^o and Ao- The derived 
values except for A478 are consistent with h = 0.7 ± 0.1 from the HST Cepheid distance 
(Tanvir et al. 1995). Note that if one adopts incorrect cosmological models, the result 
would show an artificial z-dependence of H (Suto, Suginohara & Inagaki 1995). Since 
CL0016+16 is at the highest z in our sample, we plot three different estimates of Ho using 
different sets of j3 and 9 C from radio and X-ray observations (in open pentagon, open 
triangle and filled triangle). The results are quite sensitive to the choice indicating that the 
higher resolution image analysis both in radio and X-ray observations is highly desired (see 
aslo Neumann& Bohringer 1996)0. 

Finally Figure ^| compares the angular diameter distances d A estimated from the SZ 
effect with several theoretical ones. Short dash - long dash line plots the approximation 
of d A = czif ~ 1 (l — 3z/2) valid only for z 1 and qo = 0, simply to illustrate how it 
deviates from the correct one; clearly at z £ 0.1 one should treat fl and A differently and 

3 In fact, a recent detailed analysis using the ROSAT image for CL0016+16, especially 
taking account of the ellipsoidal shape and its temperature structure, indicated a 
substantially smaller value for H than we derived here (Birkinshaw 1996). 
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the constraints cannot be given only in terms of q . Also plotted (thin solid line) is the 
Dyer-Roeder (or empty beam) distance in the fi = 1 and A = model (Dyer & Roeder 
1972) so as to indicate the possible effect of the inhomogeneity in the light propagation 
(compare with the filled beam case in thick solid line); in the light of the current systematic 
and statistical errors of the SZ measurement, the latter effect is safely neglected. 

4. Discussion and Conclusions 

We have estimated the Hubble constant H from seven SZ clusters by taking proper 
account of the effect of the cosmological parameters, which has been neglected previously. 
The estimated H varies around 30 percent depending on the assumed Qo and Ao for the 
two highest-^; clusters. Considering the uncertainties of the SZ effect (Birkinshaw et al. 
1991; Inagaki et al. 1995) and the observed scatters in the data, however, it is definitely 
premature to draw any strong conclusion at this point. Nevertheless it is reassuring that 
our estimates are consistent with Cepheid-based estimate h = 0.7 ± 0.1 (Tanvir et al. 1995) 
within 2a error-bars. In addition, the d A vs. z diagram from the future SZ measurements 
may provide a promising insight into the values of cosmological parameters as the Hubble 
diagram of distant SNe la (Perlmutter et al. 1996). 

We thank Mark Birkinshaw for various useful comments on the earlier manuscript. We 
are also grareful to Steve Allen, Akihiro Furuzawa, Mike Jones, Ian McHardy and Steven 
T. Myers for correspondences on the X-ray and radio measurements of the clusters that 
we analysed. This research was supported in part by the Grants-in-Aid by the Ministry 
of Education, Science, Sports and Culture of Japan (07CE2002) to RESCEU (Research 
Center for the Early Universe). 
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Table 1: Data for the SZ clusters 
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*) our adopted values. Our Hq is estimated for the Oo = 1 and Ao = model. 1) Yamashita 
(1994), 2) Birkinshaw et al. (1981), 3) Carlstrom, Joy, & Grego (1996), 4) Rephaeli (1995), 5) 
Elbaz, Arnaud & Bohringer (1995), 6) Arnaud et al. (1992), Wilbanks et al. (1994), 7) Birkinshaw 
& Hughes (1994), 8) Myers et al. (1996). 

a) the observed radio frequency v = 28.7GHz, b) v = 136GHz, c) v = 20.3GHz, d) v = 32GHz, 
e) in this cluster we use S* x (0) = (5.2 ± 0.83) x 10~ 6 ergcm~ 2 s -1 arcmin^ 2 directly in eq.(9) to 
compute Hq and d A , f) estimated assuming = 1 and Aq = 0. 
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Fig. 1. — H estimated from the SZ effect as a function of f2 in the case of A = (left 

panels) and Ao = 1 — (right panels) assuming the parameters in Table 1. Upper panels 
are for CL0016+16, while lower panels for A2163. The solid and dotted lines represent the 
central value and ±lcr error ranges, respectively. 
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Fig. 2.— H estimated from the SZ effect for seven clusters with la error bars in 
models with (fi , A )= (1.0, 0.0), (1.0, 1.0), (0.0, 0.0) and (0.0, 1.0). For CL0016+16, 
we plot values with different (5 and 9 C : f3 x = j3 r = 0.5, 9 X = 9 r = 0.5' (open pentagon), 
p x = 0.5, /3 r = 2/3, 9 X = 0.5', 9 r = 0.7' (open triangle), and (3 X = (3 r = 2/3, 9 X = 9 r = 0.7' 
(filled triangle). The error bars for the latter two are artficially tilted to avoid confusion. 




Fig. 3. — The angular diameter distance d A estimated from the SZ effect as a function of z 
for seven clusters (symbols). Lines indicate theoretical curves for several sets of cosmological 
parameters (h, fi , A ). The meanings of the symbols are the same as Fig. 2. 



